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{S) Abstract 

q Motivated by our research on pedestrian flows, we study a non-conservative measure- 

valued evolution problem posed in a finite interval and explore the possibility of imposing 
a flux boundary condition. The main steps of our work include the analysis of a suitably 
scaled regularized problem possessing a boundary layer that accumulates mass and de- 
tailed investigations of the boundary layer by means of semigroup techniques in spaces of 

i— i measures. We consider passage to the singular limit where thikness of the layer vanishes 

(resembling the fast reaction asymptotics typical for systems with slow transport and 

P_h rapid reactions). We obtain not only suitable solutions to the measured- value evolution 

(-h problem, but also derive a convergence rate for the approximation procedure as well as 

the structure of (flux) boundary conditions for the limit problem. 
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1 — 1 1 Introduction 

This paper discusses the resolution of various analytical issues that arrise when one wants 
to introduce boundary conditions in a measure-valued formulation of an interacting par- 
t-H ticle system or population dynamical model, in particular flux boundary conditions. The 

non-smoothness of the setting makes it unclear how to formulate such conditions, since 
the concept of normal derivate cannot be used readily. Even more important from appli- 
cation point of view - which we shall discuss below - is to understand how the boundary 
conditions, once obtained in any meaure- valued formulation, derive from particular types 
of interaction of particles or individuals at microscopic scale with the boundary in a thin 
• • boundary layer. This derivation will involve a limit in which the thickness of this bound- 

m ^ ary layer shrinks to zero. A natural problem to consider then is to determine in what 

sense and how well the limit system approximates the microscopic model with boundary 
5_h interaction layer. 



Searching for the correct flux boundary conditions is a topic often addressed in the 
literature, but this has not yet been treated in the measure-theoretical framework. We 
refer here for instance to [19j [271 [28] (in the context of reaction and diffusion scenarios) 
and Gurtin |21j (the shrinking pillbox principle in continuum mechanics) . Our motivation 
stems from insisting on capturing both the direct effect and feedback of granularity (i.e. 
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presence of a finite set of interacting individuals in a background medium) or of some 
local fluctuations (local perturbations by random forces) on the meso- and macro-scale 
dynamics of pedestrian flows. As practical goal, we search for the mechanisms responsible 
for producing streaming flow patterns able to fluidize locally congested pedestrian flows. 
Transporting measures is a tool that perfectly matches this goal especially if one has in 
mind infinite domains |13j . What happens at and close to active boundaries (e.g. where 
Cauchy fluxes need to be defined) is of primary interest in building design and evacuation 
control, but mathematically, coping with boundary (flux) measures in a system of balance 
laws is very much open. 

Measure- valued formulations gained interest substantially lately, both from fundamen- 
tal mathematical perspective (like gradient flows in metric spaces @]) as in applications 
(like evolution of structured cell densities [TJ El HU [TUl [22] > crowd dynamics 29, 6 ). Pedes- 
trian crowds and their dynamics in high-density regimes is a modern topic of intense study 
not only in security, logistics, and civil engineering (crowd-structure interactions) but also 
in non-equilibrium statistical mechanics of social systems; see e.g. |23j (an example of a 
very influential paper in the field), and [31] for a detailed overview of the current status 
or research, and applied mathematics cf. e.g. G21 [HJ [8] and references cited therein. 

In order to focus on the essential mathematical issues, we consider in this paper the 
simplest system that we can envision: a system of particles moving in a one dimensional 
confined space, say with position x in the interval [0, 1], forced by an externally determined 
velocity field v{x). There is no interaction (yet) among individuals and the boundaries are 
' sticking 1 and partially absorbing: once a particle arrives at the boundary its stays there 
and can be removed from the system (being 'absorbed' or 'gated') randomly at a time 
after arrival that is exponentially distributed with constant absorption rate a. Formally, 
one would expect an equation of the form 

d & 

Q i »t + ^ c {v»t) = -a l i t ({l})6i. (1.1) 

That is, in the measure-valued formulation for the associated particle system with 'stick- 
ing' boundary conditions [33 these conditions should be incorporated in the measure- 



valued equation (1.1) as a density dependent (point located) sink. 

There is an issue with the precise interpretation of ( fTTTj ) however. Evaluation of the 
measure at the set {1} seems to inhibit a proper interpretation in distributional sense 
on E + x [0, 1], which is the point of view taken in e.g. [TO]. The main source of 



problems in the analysis of ( 1.1 ) is the discontinuity of the map /i h- ^ /i({l}) in the natural 
topology on measures for this problem. 

Our approach differs from those cited above. It relies on our related preliminary studies 
[18111 1| (crowd dynamics, balanced mass measures) and |26| (measures in the semigroups 



context). Concerning the interpretation of (1.1), we view it as formal expression that the 



solution is obtained as perturbation of the semigroup (Pt)t^o of mass transport operators 
on finite Borel measures on [0, 1] along characterists defined by the velocity field v by 
means of a perturbation /i i— > — a/x({ 1})<5i. The solution should satisfy an integral equation 



that may be considered as a variation of constants formula for (1.1) when viewed as an 
evolutionary equation: 



Pttio - a J ds • 8l (1.2) 

o 

(see Section [2] for further discussion). 

Our approach connects to the approach to structured population models as described 
for instance by Diekmann & Getto [H] , Gwiazda et al. [221 02] or Ackleh et al. [TJ [2] • 
See also Canizo et al. [5] that simplifies part of the approach in [531 EH] ■ Their approach 
consists of first studying the situation of independently moving individuals described 
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by a (semi-)linear system like Then use the obtained results to treat the case of 

forced velocity fields changing in time and finally close the system for a density dependent 
velocity field, which models interacting individuals. This paper addresses the first step in 
this 'program' for a system with boundary conditions. Work on the remaining steps is in 
progress. 

At points in this paper the reader may find some proofs overly complicated for the 
particular one-dimensional example that is discussed. However, where possible we have 
chosen on purpose to give arguments that employ the general structure of the problem 
rather then the particular. These proofs may then be used as 'template' for dealing with 
the more delicate higher dimensional cases (2D for crowd dynamics, 'any' for structured 
population models) or abstract systems on general Polish spaces with non-empty bound- 
ary. 

Various technical complications will arrise in the higher dimensional setting. The 
boundary would need to be sufficiently smooth so that a boundary layer can be introduced 
as a suitable Lipschitz deformation of the boundary. The transport part is expected 
to become anisotropic, in the sense that some directions will matter more than others. 
Another generalisation of the currently discussed problem would be a system like 

dtih + div(«A*t) = -*R(Mt({l}))Q(i/t({l}))*i 

dtv t = +aR( f i t ({l}))Q(is t ({l}))6 1 (1.3) 

and corresponding regularizations, where the partial reaction rates Q(-) and R(-) have 



suitable algebraic structures. (1.3) brings us very close to a more fundamental chemical 



reaction theme - handling nonlinear chemical kinetics based on the mass-action law. The 



type of nonlinear structure in (1.3) is the main point of the setting addressed by Fibich 
et al. in [19]. 

1.1 Organization of the paper 



For convenience of the reader Section[L2]collects a series of properties and results on finite 
Borel measures, topologies on spaces of such measures and associated metrising norms 
and (Bochner) integrals involving measure- valued kernels that will be used throughout the 
paper. Section [2] sets-up the detailed mathematical model that will be considered. Section 



2.1 derives fundamental estimates for the movement of single individuals in the domain, i.e. 
the individualistic flow, where there is no absorption yet, but sticking boundaries only. 
Section |2.2| provides a probabilistic underpinning of the fundamental integral equation 



(1.2). In Section 2.3 we discuss an alternative approach that consists of extension of the 
domain to an unbounded setting that may seem a way to avoid some technicalities. We 
show that it suffers from the same type of technical difficulties. Section[3]discusses systems 
in which there is interaction in a (thin) layer near the boundary. There, well-posedness 
and positivity for measure- valued solutions is proven for these systems. Finally, Section [4] 



deals with the integral equation (1.2). Well-posedness is shown in Section 4.1 exploiting 



the convenient structure of the system. The most important parts are Section |4.2| and 



4.3 In the first the unique solution to ( 1.2 ) is shown to arrise as limit of the solutions of a 
sequence of thin boundary layer solutions, in the limit of vanishing thickness of the layer. 
Finally, Section [4. 3| computes the rate of convergence of these boundary layer solutions to 



the solutions of the limit equation. 



1.2 Preliminaries on measures 



If S is a topological space, we denote by M(S) the space of finite Borel measures on S and 
M + (S) the convex cone of positive measures included in it. This cone defines a partial 
ordering on M(S): \i < v iff v — /i e M + (S). Clearly, fi ^ v if and only if fJ,(E) < v{E) 
for all Borel measurable E C S. 
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For x G S, S x denotes the Dirac measure at x. Put 

(H,<t>) := I tf>dn, (1.4) 



the natural pairing between measures fi G AA(S) and bounded measurable functions 4>. If 
4> : — > S is Borel measureable, the pws/i forward or image measure of /i under $ is the 
measure 

:=/io 4> _1 . 

It is easily verified that </>) = o $). 

We denote by Cb(S) the Banach space of bounded continuous functions on S 1 equipped 
with the supremum norm || ■ Hoc. The total variation norm \\ ■ ||tv on A4(S) is given by 

|M|tv:=su P {m \ <f>€C b {S), ||0||oo ^ l} - (1.5) 

It follows immediately that for 4> : S — > S continuous, ||4>#^||tv ^ 1 1 A 4 1 1 t v - 

The total variation norm is too strong for our application, since \\S X — ^j,||tv = 2 if 
x 7^ y. The natural topology to consider is the weak topology induced by Cb(S) through 



the pairing (1.4). In this topology x h-> S x : S — > M + (S) is continuous. 

In our setting, S = [0,1] is a Polish space. It is well-established (cf. [Ml [17]) that 
in this case the weak topology on the positive cone A4 + (S) is metrisable by a metric 
derived from a norm, e.g. the Fortet-Mourier norm or the Dudley norm, which is also 
called the dual bounded Lipschitz norm, that we shall introduce now. To that end, let d 
be a metric on S that metrises the topology, such that (S, d) is separable and complete. 
Let BL(S', d) = BL(S') be the vector space of real-valued bounded Lipschitz functions on 
(S,d). For 4> G BL(S), let 

l0k:=sup{M|^ K , eWy } (1 . 6) 



be its Lipschitz constant. Then 



BL 



\<P\l (1.7) 



defines a norm on BL(5) for which it is a Banach space. In fact it makes BL(S') a Banach 
algebra for pointwise product of functions: 

||<^||bl< IMIblMbl- (1.8) 

Let || • ||b L be the dual norm on the dual space BL(5)*, i.e. for any x* g BL(5)*: 

||xl BL := sup {| (**,</>) | UeBL(S), Wbl<1}. (1-9) 

The map fi 1— > 1^ with 1^(<^>) := (fJ-,4 1 ) defines a linear embedding of M(S) into BL(S*)* 
(|16|. Lemma 6). Thus || • || BL induces a norm on A4(S), which is denoted by the same 
symbols. Generally, 1 1 A* 1 1 bl. ^ 1 1 A*l I tv - For positive measures the norms coincide: 

||/i|| BL = n(S) = HHItv for all /i G M+(S). (1.10) 

The || • || BL -norm topology on Ai + (S) coincides with the restriction of the weak topol- 
ogy ([IB]. Theorem 12). M(S) is not complete for || • || BL generally. We denote by M(S) B l 
its completion, viewed as closure of A4(S) within BL(S')*. A4 + (S) is complete for || • ||bl, 
hence closed in M(S) and M(S)bl- Note that BL(5, d) will vary with d, hence || • || B l 
on A4(S) will depend on d too. The subscript 'BL' in the notation for the completion 
^4(S')bl is intended for reminding to this dependence. 
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The || • ||g L -norm is convenient also for integration. If (X,H) is a measurable space, 
then a map x i— > /(x) : X — > M + (S)bl is Bochncr measurable (as a function mapping 
into the Banach space M(S)bl) iff a; h- f(x)(E) is measurable for each Borel set E C S 
(cf. [25] , Proposition 2.5). If z/ is a finite positive measure on (X, £), x ^ f(x) is Bochner 
measurable and x t— > ||/(x)||g L is integrable with respect to v, then / is Bochner integrable 
and 



f(x) dv(x)\ 



BL 



f(x)\\ Bh du(x) 



(1.11) 



(see e.g. [IS])- Moreover, for any continuous map P : M + (S)bl —> M + (S)bl that is 
additive and positively homogeneous, i.e. P(afi) = aP(fj,) for a ^ 0, one has 



P(j f(x)dv{xj) = J P[f{x)]dv{x), 
s s 

f{x)dv{x)){E) = J f{x){E)dv{x) 
s 



(1.12) 
(1.13) 



(cf. [25], Proposition 2.6 for (1.13)). Note that the continuity of the map x H > 5 X : S 
A4 + (S)bl together with (1.13) yields the useful identity 



A* = / $x dfj,(x). 



(1.14) 



2 Model formulation 



An 'educated guess' may lead directly to equation (1.1). In cases where the dynamics is 
more complicated - in higher dimensions, in the interior of the domain, at the boundary - 
it may not be that evident. In this section we explain how ( 1.1 ) can be derived by consid- 



ering a thin boundary layer in which well-described interactions with the boundary occur, 
and letting its thickness vanish in a limit. Such an approach would allow to determine 
the manner in which boundary interactions should be incorporated in the measure- valued 
formulation. In [TH] (in a non-measure setting of concentration functions and surface den- 
sities) such an approach was used to establish the right boundary conditions for chemical 
reactions occuring on a surface in a three dimensional domain. 

First, the operator fi i— > —-^(vfi) may be considered as the generator of a strongly 
continuous semigroup in A^([0, 1])bl: the semigroup (Pt)t^o of mass transport along char- 
acteristics associated to the velocity field v, which we shall define rigorously in Section [2~T| 
Then (1.1), when considered as an evolutionary equation, can be viewed as perturbation 
of (Pt)t^o by means of F : ^ >-» — a/i({l})<5i. The resulting perturbed solution should 
satisfy a variation of constants formula of the form 



P t ^ s F(p s )ds. 



(2.1) 



Since 5i is invariant under (Pt)t>0) i-e. Pt&i — Si for all t, equation (2.1 ) reduces to (|1.2 



Thus, as (mild) solution to ( |1.1[ ) we consider any continuous map fi : R + — > M.([0, 1])bl 
that satisfies ( 1.2 ). Note that the integral in ( 1.2 ) or (2.1 ) is well-defined, since t H> /it({l}) 
is measurable if t i-> fit is continuous for the || • ||g L -topology (see Section 1.2 1. These 



integral equations will be the main object of study in the remaining part of the paper. 



The main issue with the perturbation formulation ( |2.1| is, that the map fi t— > /i({l}) is 
not Lipschitz continuous, not even continuous on A4([0, 1])bl- The standard arguments for 
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solving such equations use a technique similar to Picard interations, but require Lipschitz 



continuity of t he 
In Section 



4.2 



perturbation term to invoke Banach Fixed Point Theorem, 
we consider F as a limit of maps that are Lipschitz continuous, 
using interactions within a boundary layer instead of solely at the boundary. The solution 
/4 to the regularized integral equations corresponding to F^ converges to the solution 



of (1.2 1 as n —> oo (see Section 4.2). Integral equation ( |1.2[ ) is closely related to a 
probabilistic description of the system, see Scetion |2.2| In Section |2.3| we show that an 
alternative approach by means of extending the state space [0, 1] to R + meets similar 



issues as (1.2| 



2.1 Mass transport along characteristics 

We assume that a single particle ('individual') is moving in the domain [0, 1] determinis- 
tically, described by the differential equation for its position x(t) at time t: 

x(i) = v(x(i)), , , 

x(0) = x , ( ' 



where v : [0, 1] — > R is a Lipschitz function. Thus, a solution to (2.2) is unique, it exists for 
time upto reaching the boundary or 1 and depends continuously on initial conditions. 
Let x{-;xq) be this solution and I Xo be its maximal interval of existence. Put 

T~d(x ) := supI Xo G [0,oo], 

i.e. the time at which the solution reaches the boundary (if it happens). 

In order to focus on the essence of mathematical issues that arrise when having a 
sticky, partially absorbing boundary, we focus on a situation where there is only one such 
boundary and this boundary point can be reached. That is, we assume v(0) > and 
v(l) > 0, such that the boundary point 1 is the only one of interest. 

The individualistic flow on [0, 1] is the family of maps $t : [0, 1] — > [0,1], t> 0, defined 

by 

^-{f^' if i €4 °' (2-3) 
I 1, otherwise. 

Lemma 2.1. (&t)t>o * s a semigroup of Lipschitz transformations of [0,1]. Moreoever, 
ft) |#tU<eH**/orOO. 
(ii) For any t, s G R + , 

sup \$ t (x) - $ s (x)\ < Hoo|t-*|. (2.4) 

xe[o,i] 

Proof, (i): The solution x(t) = x(t;Xo) is given (implicitly) as solution to the integral 
equation 

t 

x(t) = x + J v(x(s)) ds, for all t € I Xg . (2.5) 
o 

Application of Gronwall's Lemma thus yields 

|$ t (z ) - $tOr' )K \x -x' \e^ Lt for all t G I XQ n I x - . (2.6) 



A complication in proving (i) arrises from the validity of (2.6) for t in a set that depends 
on the initial conditions. Since we assume that v(l) > 0, there exists e > such that 
1 1 v $t(xo) is strictly increasing for 1 — e < x < 1. If x G (1 — e, 1) and x' = 1, then by 
monotonicity, 

|$ t (x ) - 1| < \x - 1\ (2.7) 
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for all t ^ 0. Note that there exists r > such that for any xq , x' £ [0, 1 — e] , [0, r] c I Xo H 
Jx^. Thus (2.6) holds for all < t < r and xo,x' a €[0,1 — e]. The case xo,x' £ (1 — e, 1) 
remains. We may assume Xo < Xq, hence Zj/ C Z^. Put t := su Pxo an d := supZ^. 
Then t' < to- For ^ t ^ t , estimate (2.6) holds. For t' < t < to, one has according to 
dl7l 



iStOco) ~ ^K>)l = l^-t^M - IK |** (*o) - *f (4)l 
< |x -x |eH^o <j |^o — ^|el^l-* 



For t ^ t , $t(zo) = ^t(x' ). Thus we conclude that (2.6) holds for all t £ [0, r] and 
Xo,x' Q £ [0, 1]. By the semigroup property 3>t+ s = 3?t ° each <I> t is Lipschitz and for 
n G N large such that t/n < r, 

|*tk < |*t/»|2 < eM**. 
(m): Let t, s G Z Xo . Without loss of generality, assume that t > s. 



\x(t) -x(s)\ = | J v{x(a)) da - J v{x(a)) da\ j\v(x(a))\da 


< IMU (t - S). 



(2.8) 



If both i, s £ M. + are not in I XQ , then inequality (2.8) is trivially satisfied. Suppose now 
that s £ I xo , while t is not. Let to := supZ^. Then 

\x(t) - x(s)\ = |1 - x(s)\ = \x{t ) - x(s)\ s? \\vWov (t - s), 



according to (2.8). Clearly t — s < t — s. The estimates are independent of x € [0, !]■ 
Thus we obtain fl2.4}. □ 



Define P* to be the lift of $i to the space of finite Borel measures 1]) by means 

of push forward under <£> t . That is, for all fi £ M([0, 1]), 



(2.9) 



Clearly, P t maps positive measures to positive measures and Pt is mass preserving on 
positive measures. Since the maps t ^ 0, form a semigroup, so do the maps Pt in the 
space .M([0, 1]). That is, (Pt)t^o is a Markov semigroup on [0,1]. One has ||P t /i||xv ^ 
1 1 A* 1 1 tv f° r general fx £ M([0, 1]). 

Lemma 2.2. Let fj, £ M([0, 1]) andt,s£R+. Then 

(i) UPm-PsmIIbl < IHIoo|NlTv|*-a|. 

(ii) \\P t (i\\* BL < max(l, |* t | L ) IImII^ < e^M BL . 

Proof. For all <j> £ BL([0, 1]), 

| (Ptll- P s jJL,4>) I = I 0*, O * t — O *,) | < HmIItV 11^ ° $t - 0° * S ||oo 

< IImIItv Mz su p l$t(z) - ^0*01 < IImIItv Hz IMU I* — «1, 

se[o,i] 



where we used Lemma 2.1 (ii) in the last inequality. For the operator norm of Pt on 

-\4>\l\Hl)- 



M([0, 1])bl, use that for any <j) £ BL([0, 1]), 

\{p t ii,ct>)\ = \{ii,4>o<s> t )\ <: < IImIIbl(II0I 

The statement in the lemma follows. 



□ 
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We shall need the following result on the time of arrival at the boundary, tq{x). Put 
Sq := {x £ [0, 1] | tq(x) < oo}, the set of all points that will reach the boundary in finite 
time. 

Lemma 2.3. If v is Lipschitz and v(0) > and v(l) > 0, then Sq is non-empty and 
open in [0, 1] and for each x £ Sq there exists and open neighbourhood U of x such that 
T d\u '■ U — > K is Lipschitz. 

Proof. Sq is the connected component of the set {x € [0,1] \;v(x) ^ 0} that contains 
1. So it is non-empty and open. Let v be a Lipschitz extension of v to K + such that 
\v\l = Ml and ||w||oo = IM|oo- Denote the globally existing solutions to x'(t) — v(x(t)) 
with initial value x by x(t;x). The flow map (t,x) i— > x(t;x) is locally Lipschitz, because 
v is Lipschitz. tq(x) is implicitly defined by x(tq(x); x) = 1 for x £ Sq. Since v(l) ^ 0, 
an implicit function theorem for Lipschitzian maps (e.g. [T^], Section 7.1 or [3D]) yields 
the type of local Lipschitz continuity of tq on Sq as stated. □ 

We will be interested in the part of the time interval [0, t] that an individual spends in 
the interior of the [0, 1] before possibly sticking to the boundary at 1. This time is given 
by 

tq(x) At:= min(T d (x),t), 
when the individual started at x at t = 0. 



Corollary 2.4. Assume the conditions of Lemma 2.3 Then x i— > tq(x) At is Lipschitz. 



Moreover, 1 1— > |r^(-) A i|i is non- decreasing. 1 1— > tq(x) At is Lipschitz for fixed x £ Sq. 

Proof. The set if t := e [0,1] |;rg(x) ^ t} is compact, hence it can be covered by 
finitely many open sets U as in Lemma |2.3[ So tq is Lipschitz on this set. The function 
tq{x) A t is a Lipschitz extension of this restriction to the full set [0,1]. The sets K t 
are increasing with increasing t, therefore the Lipschitz constants are non-decreasing. For 
fixed x, 1 1 — ^ tq(x) At is the infimum of two Lipschitz functions tq(x)\ and t (as a function 
of t), hence is Lipschitz. □ 

2.2 A probabilistic interpretation yielding the integral equation 

Let us consider in this section N individuals in a confinded space, with position X\ G 
[0, 1] at time t say (i = 1, . . . , N). We assume there is a sticky boundary at 1 that is 
partially absorbing. By this we mean that at the absorbing boundary we have a 'gate' 
that absorbs an individual present there a time T after arrival, which is an exponentially 
distributed random variable with (constant) rate a. We assume that the individuals are 
indistinguishable and the absorption of individuals (gating) occurs independently. We 

(i) 

denote by 7r t the law of X\ when Xq is distributed according to the probability measure 

7T . 

Since individuals behave independently from each other, the expected number of in- 
dividuals in a Borel set E C [0, 1] is given by the measure ^t(E) , where Ht satisfies 



Ht(E)=E 



The measures nj. should satisfy 



N 



N 



$>rw (2.io) 



7r t (E) = PtME) - S^E) [ an s ({1}) ds ■ <5 1; (2.11) 



s 



such that (2.101 together with (2.11) yields (1.2). To see (2.11), let us introduce the 
conditional probability 



p(t, At) := Prob (individual is gated in [t, t + At] \ X t = l) 



1 - e 



-a At 



(2.12) 



Discretize the time interval [0, t) into v steps of length As{ and let Sj be the left point of 
the i-th subinterval. Then the probability that the individual has been gated in [0, t] is 
approximately 

5>( fli ,A«K({l}) = ^ ? ' (s ^ ) ^({l})A Si 

i=i i=i 

t 

— > J air s ({!}) ds as v — > oo. (2-13) 



Formulating now (2.11 ) in terms of measures living in .M + ([0, 1])bl gives 

t 



71* = -PtTTO - y OTTs ({!}) ds. 



(2.14) 



2.3 An alternative approach by extension 

It may seem feasible at first sight to use the extension procedure as used in the proof of 



Lemma 2.3 That is, extend the velocity field v on [0, 1] to a field v on K + such that 
IMloo = IMloo and \v\l = \v\l- The solutions x(t) = x(t; xq) to 

x'(t) = v(x(t)), x{0)=x o eR + 

exist for all time and define a dynamical system (<tt)t^o consisting of the solution maps 
$t(xo) := x(t; x ). Lift $ t to measures by means of push forward: P t fj, := $t#/z. We may 
then project all excess mass in [1, oo) onto 1 by means of push forward under the map 
# : R+ -> R+ : x i-J- max(x, 1): 

Pt ■= *#A 

Pt does not depend on the choice of extension v, because v(l) > 0. It coincides with the 



mass transport operators along characteristics of Section 2.1 

One could think, that the extension procedure would also work for describing the 
partial absorption at the boundary. That is, suppose that at each point x ^ 1 there is 
absorption of mass at exponential rate a. A solution to the extended system, p, t , then 
should satisfy the variation of constants formula 



A** 



Pt^o - J P-s [al[i,oo) • fis] ds, (2-15) 



where the integral should be interpreted as Bochner integral in A^(S')bl taking values in 
Ai + (S). Since "J is Lipschitz, the map fi i-> *£#n is continuous on M + (S)bl- Thus, by 



applying the latter map to both sides in (2.15) and using ty#P t [i = P t (ty#fi), one obtains 



- Pm - / P t -.[o(*#/i.)({l}) ■ tfi] ds. 
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That is, m := x J , #^t satisfies (2.1) 



Apparently, one may consider (2.15) instead of (2.14) with the mentioned technical 



complications. Integral equation (2.15) suffers from similar problems however. The per- 
turbation F:|i4 — alji.oo) ■ [i is continuous as map A^S^tv — ► A4(5)tVj but it is not as 
map M(S)bl —> A / J(S')bl- This is easily seen by considering a sequence x n f 1: F5 Xn = 
for all n, while F8\ = — a8\. A serious issue is that the maps t i-> P t fi : R + —> A4(S)tv 
are not Bochner measureable for all /i € Ai(S). If they were, then these maps would 
be continuous on (0,oo) (cf. Theorem 10.2.3, p. 305). This cannot hold for Dirac 

measures. So ( 2.15[) cannot be interpreted as an integral equation in A4(S)tv- It can 
as an equation in A4(S')bl, but then failure of (Lipschitz) continuity of the perturbation 
leads to similar issues as for (2.14) in establishing solutions using 'standard' techniques 



for solving such integral equations. 



We prefer to consider solutions to (1.2) instead of (2.15). In the latter situation, one 



needs to consider a somewhat artifical extension to a larger state space. In fact, in case 
of a dynamical system (<&t)t^o is a general Polish space S, there is no canonical candidate 
for the extension space in which S embeds. It seems more natural then to use arguments 
'within' S. 



3 Regularized systems 



The main problem with integral equation ( 1.2 ) (or ( 2.15 )) is caused by failure of continuity 



of the map restricting measures to a measurable subset E: fi i— > 1^ • /i. Multiplication by 
a bounded Lipschitz function instead of the discontinuous indicator function improves its 
properties for the || • ||g L -norm topology: 



[0,1] 



be measureable and define Afy(/i) := f d\i for \l £ 



Lemma 3.1. Let f 

M([0,1}). 

(i) If f is bounded, then (/i)|| xv s$ H/IUII^IItv- 
(ii) 7//eBL([0,l]) ; then \\M f (fi)\\* BL < ||/||bl||mIIbl- 
Proof. Part (i) is easily verified. For part (ii), let <f> € BL([0, 1]). Then 

\(M f {^),4>)\ = | <m, / • 0> | < IImIIblII/ • IMIblII/IIb L ||< 

because BL([0, 1]) is a Banach algebra. 



BL: 



□ 



by 



The integral equation (2.1) in which the perturbation F : fj, <-> afi({l})Si is replaced 



Ff : fi i-> af(-) djj, 



with / G BL([0, 1]) such that ^ / ^ 1 and /(l) = 1 will be called a regularization of 



( 1.2 1. In a regularized system particles can be absorbed also in the interior of the interval 



where / > 0. The function / is called the regularizer. 



3.1 Well-posedness and positivity of solutions 

We consider integral equation ( |2.1[ ) with F = Mf for / e BL([0, 1]). 

Theorem 3.2. Let p, e M([0, 1]) and f e BL([0, 1]). Then there exists T e (0, oo] and 
a unique solution [i t to integral equation (2.1) with F — Mf in C([0, T), Ai([0, 1]bl)- T 
is finite if and only if \\fi{ ||xv — > oo as t — > oo. The solution depends continuously on 
initial data. 
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Proof. In [9], Section 2.2, the case of well-posedness for a non-linear and time-dependent 
perturbation of a transport equation on M. d in the space of measures is considered. The 
line of proof in our bounded setting is completely similar, hence is omitted. Note that 
our pertubation F^ n > is linear and bounded for || • ||tv, such that it satisfies the crucial 
condition (H5) in [9] that ensures solutions starting in /zo € M.([0, 1]) will remain measures 
up to time T where blow-up occurs in || • ||tv (|2]i Theorem 2.4). □ 

Theorem 3.3. Assume that fi e A^ + ([0, 1]) and let f £ BL([0, 1]). The solution fi{ to 
(2.1 1 with F — Mf exists for all time and /if € .M + ([0, 1]) for all t. 



Proof. We will show below that the solution remains a positive measure on its domain of 
existence [0, T). Then we can apply Gronwall's Inequality, since for all t g [0,T), 

t 

II^IItvHI/^IIbL<II^oIIbL + / \\Pt- S M f (^)\\* Bh ds 



f 

<eH-*|| Mo || TV + J e M^-s)||/|| BL || M {|| TV d S . 
o 

Therefore 

||m{||tv *S ||Mo||Tve (Mi + ll/l|BL)t 

and the total variation norm cannot blow-up in finite time. So T = oo, according to 
Theorem E21 

(Positivity). The arguments used in [9], Lemma 3.2, to prove positivety cannot be trans- 
fered simply to our setting, since they use reduction to i 1 -solutions and approximate 
measure-valued solutions by the former. However, in the bounded domain there are no 
-^-solutions (for all time). The idea of proof is essentially the same however. Let / be 
the positive and negative part of /, such that / = / + — /~. Put b := ||/~||bl- The per- 
turbation Mf + bl is then a positive perturbation, mapping .M + ([0, 1]) into itself. Let A 
be the generator A of P t on A^([0, 1])bl- Then A — bl is the generator of a C -semigroup 
too: that of the positive semigroup Pt := e~ bt P t . One can show (e.g. [33], in quite a 
general context) that a solution to the integral equation 



^ = P tMo + / Pt-s(M f + bl)(» s ) ds (3.1) 



is a solution to (2.1) with F = Mf (and vice versa). The perturbation in (3.1) is positive 



and Lipschitz for || • ||g L . Hence (3.1) can be solved (on a suitably small time interval 



[0,Ti]) by means of Picard interations: 

u (t) := P t n , u k+ i(t) := P t ^ + I P t - s {Mf + bl)(u k (s)) ds. 



o 



Since uo(t) is a positive measure for each t, each Uk{t) is, and positivity follows on [0, n]. 
This solution can then be extended to [ti,T2] as positive solution. Thus continuing one 
obtains a positive solution on [0,T). □ 



According to |3.3| for / 6 BL([0, 1]) we can define a semigroup of solution operators 
?f/x :=/i< onM+([0,l]). 



11 



Lemma 3.4. The solution semigroup consist of Lipschitz maps on A^ + ([0, 1])bl' 

HQ^o - Q{moIIbl < e< Hi+ll/l|BL >||Mo - A*oIIbl- 

Proof. The proof consists essentially of repeating the type of estimates and arguments 
employed in the proof of the first part of Theorem 3.3 so these will be omitted. □ 

3.2 Boundary layer approximation by regularized systems 

We shall now consider a countable family of regularized systems defined by a decreasing 
sequence (/„) C BL([0, 1]) of regularizers. Define f n {x) ■— [n(x — (1 — ^))] + , where [■] + 
denotes the positive part of the argument. That is, 

f(x)-l°> for xe [0,1 -^); (32) 



n(z-(l-i)), /orie[l-il] 
It is easy to see that 

||/n||BL = ||/n||oc + |/nU = l + n. (3.3) 

Moreover, 

Lemma 3.5. For any m, n G N + , satisfying m ^ n 

\\fn ~ fmWoo = 1 - (3-4) 

m 

Define : M([0, 1]) -> M([0, 1]) by means of 

FW/i = -a/ n -d/*. (3.5) 

Lemma 



3.1 



and (3.3) yield that is a continuous linear operator on Ai([0, 1])bl with 
operator norm \\F&>\\ < |a|(n + l). The norms of the F<") are not bounded in n unfortu- 
nately. This complicates the transition to the limit equation as n — > oo and necessitates 
the more delicate approach presented in Section [4] below. However, if measures are prop- 
erly ordered, then one has 

Lemma 3.6. For any n G N + and /x, v G A / J + ([0, 1]) satisfying fi ^ v, the following 
inequality holds: 

||i^-F (n M B L<a||/x-H| B L- (3-6) 

Proof. By assumption, fi - v G M + ([0,1]), so -fW^ - i/) e 7V4 + ([0, 1]). Therefore 
— f ||bl = Wl 1 ~ HItv and Lemma 3.1 (i) yields 

WFWp - F™u\\*bl = \\F {n) (v ~ v)\\tv < - HItv = a\\p - v\\* Bh . 

□ 

We shall denote by /4 the mild solution to the regularized system defined by /„ 
with initial condition /l*o £ A4 + ([0, 1]). That is, i i-> /4™' is the unique continuous map 
K+ -> A4+([0, 1])bl that satisfies 



-PtMo + y P t - s F (n Vi n) ds, for all f ^ 0. (3.7) 



(n) 

IH 



As a consequence of Theorem |3.3| and Lemma [3T4] we thus have: 

Corollary 3.7. For each n G N + and /i G A / J + ([0, 1]) there exists a unique mild solution 
G C(R + , A4 + ([0, 1])bl)) t° equation (3.7). The solution operators Qj defined by 
Qt™Vo : = A*{ are Lipschitz continuous on .A4 + ([0, 1]) /or || ■ ||g L . 
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3.3 Ordering of regularized solutions 

Lemma |3.6| indicates that ordering of the regularized solutions may be exploited in esti- 
mates. One may anticipate such an ordering, since with increasing n, the functions f n 
decrease. Thus one would expect more mass in the mild solution for larger n at the same 
time, because less absorption has taken place. We shall make this precize below. 

Lemma 3.8 ('Diracs stay Diracs'). If fJ-o = $x> then is concentrated at & t (x). 

Proof. On a small interval [0, T] the solution to /it = Pf/io + L Pt- S F^ fx s ds is given by 
the limit of a sequence u k (t) € fit € C([0, T], .M + ([0, 1])bl) defined by 

t 

«o(t) := PtMo, u*+i(t) := f*Mo + / Pi- S F (B) %(#. (3.8) 



If /io = S x , then P*/xo = 5$ t ( x j is a Dirac distribution located at $t(a;). Suppose now that 
Ukit) = Ck{t)5$ t ( x ) for some Ck(t) > 0. Then 

F (n) u k {t) = -af n ■ [c k (t)S 9t(x) ] = -af n ($ t (x))c k (t)6$ t ( x) 

is a Dirac located at $ s (x). Consequently, 

t 

u fc+ i(t) = ** t ( a .)+ J Pt- s [F (n) u k (s)]ds 



/ 

= [l-a J f n ($s{x))c k (s) ds] ■ S^ t{x y (3.9) 



Thus /i 4 := linifc-^oo u k {t) is a multiple of a Dirac measure positioned at $t(x). □ 
Substituting /j, t = a[ n <5$ t ( x ) into the variation of constants formula leads to 



Q 



Evaluating the measures on {$ t (a;)} (or on [0, 1]) yields 

t 

ai n) = 1 - a / /„($ s (x))aW( S ) ds. (3.10) 



Since i — > is continuous, it follows from (3.10) that this is even differentiable, and 
correspondingly, 

f t a^ = -af n (^(x))a^. (3.11) 



Therefore 



a[ n) = a( n) exp(- a / /„(*.(*)) ds). (3.12) 

and we obtain: 
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Lemma 3.9. A SSUUIG thdt flQ — 0Lq8 X q with ao G M + and xq £ [0, 1] . Then for any 



n £ N + , jiij 71 *' = a i n ^ , s>t(xo)' "where 



» 



= a ea;p -a / f n ($ s (x )) ds 



(3.13) 



In particular, if m,n £ N + and n ^ m, then oq ^ ctj ^ . 
This leads to 

Proposition 3.10 (Ordering of regularized solutions). For aZZ m,n € N + , suc/i i/ioi 
n m, the corresponding regularized solutions satisfy n[ n) ^ for allt^O. 

Proof. Let be the solution operator for the regularized system with regularize! - /„. 
This operator is Lipschitz continuous for || • ||g L (Corollary 3.7). Identity (1.14) and 
continuity is used to obtain 



/ xW=qW( Mo )= J Ql n \6 x )dfi (x) 

[0,1] 



tr\5 x )d^{x) = Qr i \m)^^r\ 



[0,1] 



using Lemma |3.9| 



(3.14) 

(3.15) 

□ 



Another type of ordering of regularized solutions is provided by the observation that 
measures can only loose mass over time. There is no source in the system. This yields 

Lemma 3.11. If fi e .M + ([0, 1], then for all t ) 0, n e N + , /4™' < PtfJ-o- In particular, 

ii /4 ikv < ii/^oIItv- 



Proof. If fj,Q is positive, then so is fx[ n ' (Corollary 3.7). Using the variation of constants 
formula (3.7) again, one has 



/4 n) = PtMo -a J P t _ s [/„ M i n) ] ds 



(3.16) 



The integral term in (3.16) is a positive measure, so PtHo — /4 *J 0. The norm estimate 
follows directly from the positivity of /4™' and this ordering: 



M^Htv = M r ([0, 1]) < PtlM>([0, 1]) = II^IItv = M 



,(«) 



TV- 



□ 



4 The limit equation 

In this section we consider the problem of existence, uniqueness and continuous depen- 



dence on initial conditions for the solutions of the integral equation (1.2) for the limit 
system with interaction at the boundary only. In Section |4.1| we prove this using the 
particular structure of the system that allows such an approach. In Section |4.2| consider 
the limit system truely as a limit of regularized systems with interactions in a boundary 
layer. The latter employs estimates on the integral equations for the regularized systems 
that may be used for other type of systems in a similar manner. 
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4.1 Well-posedness 

As can be expected, the results follow from careful analysis what happens to mass on 
the boundary, which is possible in our case, because of the relative simplicity of the one- 
dimensional situation. Note that the proof of continuous dependence on initial conditions 
in Proposition |4.3| requires a delicate argument, in particular it employs Corollary |2.4| 



Proposition 4.1 (Uniqueness). A solution to (1.2) in C(IR + , .M([0, 1])bl) * s unique, if 
it exists. 

Proof. In this case, a modified argument of Gronwall-type shows the uniqueness of solu- 



tions. In fact, if (1.2) had two solutions /it and (it on [0,T], having the same initial data 
//q, then for alH ^ 0, 



IM - fk = -a / [|U S ({1}) - fi s {{l})]ds ■ Si. 



(4.1) 



That is, two solutions can differ by mass concentrated at 1 only. Note that the integrand 



in (4.1) is a bounded measurable function. Evaluating the latter equation at {1} yields: 

t 

Mt({l» - £*({!}) = -a J Ml}) - Aa({l})]d*. (4-2) 
o 

and consequently 



|^({1}) -£*({!})! <a / M{l})-A.({l})|d« 



A version of GronwalPs Lemma yields that |/xt({l}) — £*({!■}) I = for all t ^ 0. 



(4.3) 



□ 



Since there is no 'smoothing' effect in the dynamics in the interior of the interval [0, 1], 
we expect Dirac masses to stay Dirac masses. These move according to Pt- The latter 
acts simply on Dirac masses: PtS x = 5$ t ( x y Therefore one may try as particular solution 
to ( 1.2 ) with /i = S x : 

IH ; = a* ( 4 - 4 ) 
Substituting (4.4) into (1.2) yields, after evaluation of the measures on the full space [0,1]: 



a x (t) = a x (0) - a j a x (s)ds, 

Tg(x)At 

where Tg(x) At is the time in the interval [0, t] that the individual is not at the boundary 
(see Section pO) . Thus we obtain that 

Ql (t)= ftl (0)e-^W Ai l. (4.5) 



Any initial measure is a superposition of Dirac masses, according to (1.14). Therefore 



we obtain the following existence result and integral representation for the solution to 



(1.21 



Proposition 4.2 (Existence). For each fiQ € -M([0, 1]) there exists a continuous solution 
pi : R+ -)• M([0, 1])bl to ([L2j) defined by 

IH'.= J e-^- T3 ^S^ {x) dfi (x). (4.6) 

[0,1] 

as Bochner integral in A / I([0, 1])bl- 
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Proof. The integrand in (4.6 ) is a bounded continuous function from [0, 1] into A4 + ([0, 1])bl 
(Corollary 2.4). Thus for /xo <E .M + ([0, 1]) the Bochner integral exists, with value in 
.M + ([0, 1]), because this cone is closed. For /i £ -M([0, 1]) the integral yields a measure in 
.M([0, 1]) C -M([0, 1])bl, by using the Jordan decomposition /i = ^ — (1q . 
Now hx t e R+ and let teR+. Then 

/ || e - o[t - T8 ^ At] (^ t(:c )-^ (x))||; L d|Mo|(a;) 

[0,1] 



\ e -a[t-rg(x)M] _ e -a[t ~Te(x)/\t ] 



d\jj, \(x) 



[0,1] 



-a[t— Tg(x)At] 



-a[t a -T g (x)At ] I 



d\no\(x) 



[0,1] 



Continuity of the maps t H ¥ <&t(x) and 1 i-» exp(— a[t — tq(x) A t]) (cf. Corollary 2.4) and 



application of Lebesgue's Dominated Convergence Theorem yield continuity of t H >• /i t . It 
is easily verified that //t satisfies (1.2). □ 



A 'standard' argument with GronwalPs Inquality to obtain continuous dependence on 
initial conditions fails in this setting, because the pertubation is not Lipschitz continuous. 



Instead we shall use (4.6) 



Proposition 4.3 (Continuous dependence on initial conditions). For each T ^ 0, there 
exists Ct > such that for all initial values /^OjMo G -M + ([0, 1]) the corresponding solu- 
tions \i and jj! to (|1.2[) satisfy 



I A** — a4IIbl < CtIImo-MoIIbl 



(4.7) 



for all t € [0,T]. 



Proof. In view of Lemma 2.2 (ii) we need to control the integral term in (1.2). Intuitively, 
this term is the total amount of mass that disappeared from the system in the time interval 
[0,t\. To be precise, according to ( |1.2[ ) and (4.6): 



i J m.({1}) ds = P t fM>(S) - fit(S) = \\fM>hv - e- at J < 



,a[Tg(x)At] 



d/j, (x) 



MoIIbl 



e- a *(^ ,e a[T8( - )A * ] 



Note that the map x h-» e a [ ra ^^ A *l is bounded Lipschitz (Corollary 2.4): 



| e a[r e (-)At]| ^ e at \T 9 {-) A tjj 



and 



a[r g (-)At] I 



s? e" 



because t (x) A t ^ t for all x. Therefore, using Lemma 2.2 and (4.8) 



(4.8) 



llM* - /4IIbl < W p t^o - Mo)IIbl + o| / /*.({!}) - M' S ({1}) d« 



^ el"l^|| M o - MoIIbl + «|HmoIIbl - IImoIIblI + e^Mo - Mo, e^^'l) 
< IIW) - MoIIbl [e [vlLt + a + e- at \\e a ^ M % L ] 
^ ( a + e l-l^ f + 1 + |T- a (-) A ||mo - MoIIbl- 



The factor in front of \\/iq — /^oIIbl * s non-decreasing in t, according to Corollary 2.4 
is dominated by a constant Ct for t in an interval [0, T]. 



so it 

□ 
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4.2 Approximation by regularized systems 



Our main point is, that solutions to ( |1.2[ ) can be viewed as a limit of solutions to regularized 
systems, with interaction with the boundary in a thin boundary layer. Well-posedness of 
the latter systems follows essentially the standard proof for abstract semi-linear equations 



in Banach spaces, see Section 3.1 Such results are expected to persist in more complex 
situations (see e.g. [9]). In this section we show how estimation of the variation of 



constants formula for the regularized systems (3.7) leads to equicontinuity of the family 
of solutions | n £ N} in C([0,T],M([Q, 1])bl), for any T ^ 0. An Arzela-Ascoli 

type of theorem [3] then yields precompactness, hence the existence of subsequences that 



converge. The limit of each convergent subsequence yields a solution to (1.2). The latter 
is unique, so we obtain convergence of the full sequence. 
First we need a lemma: 

Lemma 4.4. For all s,t £ [0,T], satisfying s ^t, the following estimates hold: 
(i) II/; Pt-vF^nP -P s - a F^^ ] da\\* BL < a|M|oo||/xo||TV«|t-»|; 
(ii) WfiPt-aFW^daW^ ^ oH/xoIItvI*-*!. 



Proof, (i): According to Lemmas 2.1 2.2 and 3.11 and properties of the Bochner integral, 
we have 

s 

|| J P t - a F^^ - P 8 _ a FV>nP da\\* BL 

Q 

s 

< J \\P-aF^^ - /V^V^IIbl da 



5 

<M|oo|*-«| /||F< n V£°llTvdff 



<o|HUI*-a| W^hv da 



^ a \\v\\oo\t - s\ ||/io||TV s - 



(ii): The second estimate is obtained as follows: 



Z I Z 

J P t ^F^^da\\ BL < J \\P t ^F^^\\ BL da = a J\\P t _ a [f n ^]\\ TW da 

s s s 

t t 
= a f\\f n ^\\ TV da^a J\\^\\ TV da 

S S 

< a IImoIItv I* - s \- 



Here we used that f n fJ,&, hence Pt-aifnf-a 1 '], is a positive measure for all a, so || ■ \\ BL 
and || • 1 1 tv coincide. □ 

We now arrive at the main result. 

Theorem 4.5. Let /iq £ A4 + ([0,l]) and T ^ 0. The solutions jif"' to the regularized 
system (3.7) converge in C([0,T], Al + ([0, 1])bl) to the unique solution /it of (1.2). 



»i 
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Proof. Let T ^ and Ho <E A4 + ([0, 1]). First we prove equicontinuity of the family of 
regularized solutions { (j, (n) \ n e N+} in C([0,T],M([0, 1])bl)- To that end, let n € N 
and t,s e [0, T] such that t > s. Then 

H^ n) -Mi n) llBL<ll^o-P s Mo||lL 

+ || J Pt^F^^da- J P s - a F^^da\\* BL 
o o 

s 

|| J Pt^FWpW - P„- a FV>nW da\\* BL 



sc HmoIItv ||f ||oo]t — s\ 



+ || J P t - s F^^da\\* Bh 

S 

HmoIItv \\v\\oo\t - s\ + IMIooallMollTV • s\t- s\ + o||/xo||tv|* - s|. 



Here we used Lemma [2. 2 1 in the first step and Lemma [4.4| in the last. Thus there exists 
C = C T > such that for all n £ N and t,s€ [0, T], 



(") 

M 



^llBL^CWxvlt-sl. 



In particular, { G N} is uniformly equicontinuous in C(\Q,T] M( \0, 1])bl)- 

For each t € [0,T] and n e N+, Hm^IItv «S ||mo||tv (Lemma 3.111. Since [0, 1] is 
compact, the set {^4"' ) | n 6 N+}, consisting of positive measures, is uniformly tight 
for each t. According to a version of Prohorov's Theorem (cf. [7J, Theorem 8.6.2), it is 
precompact in A^ + ([0, 1])bl- Ambrosetti's generalisation of the Arzela-Ascoli Theorem 
[3] yields that { | n G N+} is precompact in C([0, T],M([0, 1])bl)- 

Thus there exists a convergent subsequence in C([0,T],M([0, 1])bl) with limit 

//*. Since M+([0,1]) is closed in M([0,l]) B h, V* € C([0, T], A^+(S , ) B l)- Application of 
Lebesgue's Dominated Convergence Theorem (for Bochner integration, see [15]) to (3.7) 
shows that satisfies (2.14). Any convergent subsequence of ((jM 1 ^) must have the same 
limit, by uniqueness of solutions to (2.14) (Lemma 4.1). We conclude that the sequence 
(//(")) itself is convergent in C([0, T]^M+([0, 1]) B l) with limit fj,*. □ 



4.3 Rate of convergence of approximation 

We obtain a rate of convergence of regularized solutions to the solution of the limit 
system for a subset of initial conditions by means of a Cauchy- argument. We show that 
for ininitial conditions in this subset (which is sufficiently rich from the point of view of 
applications), the sequence of regularized solutions is actually a Cauchy sequence. 

Reconsider equation (3.7) that is satisfied by the regularized solution /4™- 1 . Let n > m. 
Then, writing || • || instead of || • ||bl*, 



IK 



(") _ ,,( m )i 



P t - S [F^nM - F< m Vi m) ] ds\ 



o 



i: / \\Pt-s[F {n) ni n) -F( n Vi m) ]ll^ + J \\P t s[F {n) ni m) -F( m Vi m) ]IMs (4-9) 
o o 



Let us denote the left and right integral terms in (4.9) by 1\ and I2 respectively and 
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start by considering the last term: 



h «S e Hi * / ||^("Vi m) - F( m Vi m) ll da, 



according to Lemma 2.2 Put E m := [1 — — , 1). 

Lemma 4.6. Lei to, n g N + suc/i t/iai n > to. TVien /or a// /x € .A4([0, 1]), 
||F<*V ~ F (m Vl|TV < a |i - -M • m| M |(£ m ). 



Proof. Let 93 € Cf,([0, 1]). Then according to Lemma 3.5 



|<F("V-F(™W)|<a y |/ n - / m | M 

<a(l-f)-||^||oo-|M|(^m)- 



(4.10) 



Note that one can remove {1} from integration in (4.10), since the integrand is zero at 1. 
The result folows immediately from (4.101. □ 

Since we take positive initial conditions, jsi" 1 ^ ^ for all s. Hence 



1 - — — 1 


•to/ 




In ml 





Now, ^ PsfJ-o (Lemma 3.11), so we obtain an estimate for I2 that depends on the 
initial condition only: 



z 

h < ae^^-^-mlJ P sl x Q ds]{E m ). 



(4.11) 



We shall now provide conditions on /zo such that the integral in (4.11 ), evaluated on E m , 
is C(1/to). Let A denote Lebesgue measure on [0, 1] 

Proposition 4.7. Let /.to € .M + ([0,1]) be such that /io = /io,s + 3o°^ * s the Lebesgue 
decomposition of fig with /io jS the part of /io that is singular with respect to X. Suppose 
go G L°°([0, 1]) and /io, s is such that for X-a.e. s 6 [0, 1], 



X{E m ) 

Then there exists c > 0, independent oft, such that 

t 



— > as to — )■ 00 . 



(4.12) 



limsup to f P sf i (E m )ds ^ cHffolUM^V^* - 1). 

m— too J 



(4.13) 



Proof. For each s € [0,t] fixed, P s fJ-o, s is singular with respect to A (cf. Lemma 3.8). 
Moreover, $ s is Lipschitz, hence differentiable almost everywhere by Rademacher's Theo- 
rem. Put J s := L>$ s e i°°([0, 1]). Then || Jsix)]]^ ^ c|$ s | L for some c> 0. Furthermore, 
for each s, 

P s {g dX) = c (s)6 1 + g s dX, (4.14) 
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where g s g ^([0, 1]) n L°°([0, 1]). Indeed, for any E C [0, 1] Borel, we have 



P s (g d\)(E) = 4> s #(g Q dX)(E) 



g dx 



go dx 



*r ({i}n£) 



c {s)S 1 (E)+ / 5o ($ s (x)) J s (x)dx, 



E\{1} 



where Co is some constant possibly depending on s. Taking g s := (go°^s) ■ J s gives (4.14|. 
Based on [20], we note that for E m := [l — i l), we obtain 



P s fj,o(Er, 



lim 

m->oo X(E m ) 



= lim 

m— >oo 

= lim 



(P fl /zo, a + co(a)(?i)(£ m ) + (g s dX)(E m ) 
X{E m ) 

P s Ho, s (E m ) (g s d\)(E m ) 
X{E m ) X(E m ) 



(4.15) 



The first term in the limit shown in (4.15 ) is dealt with by (4.12 ). Note that P s /io jS (£ , m ) < 
PsMo,s([0, 1]) < P s ^o([0, 1]). So (4.12| and Lebesgue's Dominated Convergence Theorem 
give 



X(E 



sr^r [ PsHo. s (E m )ds -> as m -> oo. 
o 



(4.16) 



Let us now turn the attention of the second term in the last limit of (4.15). For all 
s e [0,t]: 



(g s dX)(E m ) 



X{E m ) 



< \{E m ) J g ($ s (x))\J s (x)\dx 

< ||5o ° ®s\\oc,E m ■ c\$ s \ L c|| ff0 ||ooe 



So, we can conclude that 



I— I ' (g s dX)(E m )ds < cllffolloo • 7^(e Mit - 1) (4.17) 



X(E m ) 



Combining (4.16) and (4.17) together yields the desired result. 



The first integral in (4.9) may be estimated as follows: 

t 



h < e MLt / ||F("Vi n) -F("Vi m) ||ds. 



□ 



Observing that ^i™' 1 > /Lti"^ for all s e [0,i] (Theorem 3.10), Lemma 3.6 yields 



h «S ae lvlLt / \\4 n) -4 m) \\ds- 



Thus we can apply a version of Gronwall's Lemma (e.g. [35], Lemma D.2) and obtain 



1 + ate Hi *exp(ate 



(4.18) 



Here 

Thus one gets 



C t = ^||flt)||oo|«ll 1 (c Hit - 1). 
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Theorem 4.8. Let the initial datum /io satisfy the conditions of Proposition \4-% Then 
the sequence of solutions (fJ>^ n ') to the regularized systems defined by (f n ) is a Cauchy 
sequence in C([0, T], A4 + ([0, 1])bl) for each T > 0. Moreover, 

H/4 n) -Mt||! L = O(i) asrwoo, (4.19) 
uniformly on compact time intervals. 

Note t hat t he rate of convergence is 'as good as' that of the sequence (/„) in view 
of Lemma 



3.5 



A similar type of result could be obtained with regularizations /„ dif- 
ferent from f n but still bounded Lipschitz, supported on E n and such that /„ converges 
monotone (decreasing) pointwise to l{i}- The crucial condition then is, the equivalent of 



Lemma |3.5[ namely that for n > m, 

\\fn - /mlloo < u)(n,m) ■ m, 

where to ^ such that for each e > 0, there exists N £ N such that u>(n, m) < e for all 
n,m ^ N. The rate of convergence of the approximations towards the limit solution is 
then controled by uj(n) := limsup m ^ 00 cu(n, m). 
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